Piecing the barcoding puzzle of Palearctic water frogs (Pelophylax) sheds light on amphibian biogeography and global invasions.

Palearctic water frogs (genus Pelophylax) are an outstanding model in ecology and evolution, being widespread, speciose, either threatened or threatening to other species through biological invasions, and capable of siring hybrid offspring that escape the rules of sexual reproduction. Despite half a century of genetic research and hundreds of publications, the diversity, systematics and biogeography of Pelophylax still remain highly confusing, in no small part due to a lack of correspondence between studies. To provide a comprehensive overview, we gathered >13,000 sequences of barcoding genes from >1700 native and introduced localities and built multigene mitochondrial (~17 kb) and nuclear (~10 kb) phylogenies. We mapped all currently recognized taxa and their phylogeographic lineages (>40) to get a grasp on taxonomic issues, cyto-nuclear discordances, the genetic makeup of hybridogenetic hybrids, and the origins of introduced populations. Competing hypotheses for the molecular calibration were evaluated through plausibility tests, implementing a new approach relying on predictions from the anuran speciation continuum. Based on our timetree, we propose a new biogeographic paradigm for the Palearctic since the Paleogene, notably by attributing a prominent role to the dynamics of the Paratethys, a vast paleo-sea that extended over most of Europe. Furthermore, our results show that distinct marsh frog lineages from Eastern Europe, the Balkans, the Near East, and Central Asia (P. ridibundus ssp.) are naturally capable of inducing hybridogenesis with pool frogs (P. lessonae). We identified 14 alien lineages (mostly of P. ridibundus) over ~20 areas of invasions, especially in Western Europe, with genetic signatures disproportionally pointing to the Balkans and Anatolia as the regions of origins, in line with exporting records of the frog leg industry and the stocks of pet sellers. Pelophylax thus emerges as one of the most invasive amphibians worldwide, and deserves much higher conservation concern than currently given by the authorities fighting biological invasions.

Sex-biased phenotypic plasticity affects sexual dimorphism patterns under changing environmental conditions.

Sexual dimorphism is almost ubiquitous in animals. A common pattern observed across multiple taxa involves differences in development time (sexual bimaturism) and body size (sexual size dimorphism) between conspecific males and females. Furthermore, a strict association of dimorphism at these traits has been documented in several taxa, where the sex showing shorter development time also has a smaller body size than the other sex. Growth and development are strongly dependent on environmental conditions during individual life-cycle in ectotherms, inducing considerable phenotypic plasticity. However, how phenotypic plasticity affects the association between sexual dimorphism in development time and body size remains unclear. Here, we tracked development time, body size, and body mass throughout the ontogeny of the mosquito Aedes mariae. The larval development of this species is strictly linked to Mediterranean Sea rock-pools, whose highly variable environmental conditions over minimal time frames make this organism-environment system ideal for exploring plasticity-led eco-evolutionary processes. We found differential plasticity between males and females, dissolving the link between dimorphism in development time and body size under increasing temperature and decreasing salinity conditions. These findings contrast with the current hypotheses proposed to explain the origin of the association between sexual bimaturism and sexual size dimorphism, highlighting the condition dependence of sexual dimorphism patterns and the need to consider phenotypic plasticity in future studies on their evolution.

An update and ecological perspective on certain sentinel helminth endoparasites within the Mediterranean Sea.

The Mediterranean Sea is recognized as a marine biodiversity hotspot. This enclosed basin is facing several anthropogenic-driven threats, such as seawater warming, pollution, overfishing, bycatch, intense maritime transport and invasion by alien species. The present review focuses on the diversity and ecology of specific marine trophically transmitted helminth endoparasites (TTHs) of the Mediterranean ecosystems, aiming to elucidate their potential effectiveness as 'sentinels' of anthropogenic disturbances in the marine environment. The chosen TTHs comprise cestodes and nematodes sharing complex life cycles, involving organisms from coastal and marine mid/upper-trophic levels as definitive hosts. Anthropogenic disturbances directly impacting the free-living stages of the parasites and their host population demographies can significantly alter the distribution, infection levels and intraspecific genetic variability of these TTHs. Estimating these parameters in TTHs can provide valuable information to assess the stability of marine trophic food webs. Changes in the distribution of particular TTHs species can also serve as indicators of sea temperature variations in the Mediterranean Sea, as well as the bioaccumulation of pollutants. The contribution of the chosen TTHs to monitor anthropogenic-driven changes in the Mediterranean Sea, using their measurable attributes at both spatial and temporal scales, is proposed.

De novo transcriptome assembly of an Antarctic nematode for the study of thermal adaptation in marine parasites.

Understanding the genomic underpinnings of thermal adaptation is a hot topic in eco-evolutionary studies of parasites. Marine heteroxenous parasites have complex life cycles encompassing a free-living larval stage, an ectothermic intermediate host and a homeothermic definitive host, thus representing compelling systems for the study of thermal adaptation. The Antarctic anisakid Contracaecum osculatum sp. D is a marine parasite able to survive and thrive both at very cold and warm temperatures within the environment and its hosts. Here, a de novo transcriptome of C. osculatum sp. D was generated for the first time, by performing RNA-Seq experiments on a set of individuals exposed to temperatures experienced by the nematode during its life cycle. The analysis generated 425,954,724 reads, which were assembled and then annotated. The high-quality assembly was validated, achieving over 88% mapping against the transcriptome. The transcriptome of this parasite will represent a valuable genomic resource for future studies aimed at disentangling the genomic architecture of thermal tolerance and metabolic pathways related to temperature stress.

Evolution of personality and locomotory performance traits during a late Pleistocene island colonization in a tree frog.

Recent empirical and theoretical studies suggest that personality and locomotory performance traits linked to dispersal abilities are crucial components of the dispersal syndromes, and that they can evolve during range expansions and colonization processes. Island colonization is one of the best characterized processes in dispersal biogeography, and its implication in the evolution of phenotypic traits has been investigated over a wide range of temporal scales. However, the effect of island colonization on personality and performance traits of natural populations, and how these traits could drive island colonization, has been little explored. Noteworthy, no studies have addressed these processes in the context of late Pleistocene range expansions. Here, we investigated the contribution of island colonization triggered by postglacial range expansions to intraspecific variation in personality and locomotory performance traits. We compared boldness, exploration, jumping performance, and stickiness abilities in populations from 3 equidistant areas of the Tyrrhenian tree frog Hyla sarda, 2 from the main island (Corsica Island), and 1 from the recently colonized island of Elba. Individuals from Elba were significantly bolder than individuals from Corsica, as they emerged sooner from a shelter (P = 0.028), while individuals from Corsica showed markedly higher jumping and stickiness performance (both P < 0.001), resulting as more performing than those of Elba. We discuss these results in the context of the major microevolutionary processes at play during range expansion, including selection, spatial sorting, founder effects, and their possible interaction with local adaptation processes.

The ecological importance of hybridization.

Hybridization as an evolutionary process has been studied in depth over the past few decades. Research has focused on its role in shaping reproductive barriers, its adaptive value, and its genomic consequences. In contrast, our knowledge of ecological dimensions of hybridization is still in its infancy, despite hybridization being an inherently ecological interaction. Using examples from various organisms, we show that hybridization can affect and be affected by non-reproductive interactions, including predation, competition, parasitism, mutualism, commensalism, and organism-environment interactions, with significant implications for community structure and ecosystem functioning. However, since these dimensions of hybridization have mostly been revealed from studies designed to decipher other evolutionary processes, we argue that much of the eco-evolutionary importance of hybridization is yet to be discovered.

De novo transcriptome assembly and annotation for gene discovery in Salamandra salamandra at the larval stage.

Dispersal is a key process in ecology and evolutionary biology, as it shapes biodiversity patterns over space and time. Attitude to disperse is unevenly distributed among individuals within populations, and that individual personality can have pivotal roles in the shaping of this attitude. Here, we assembled and annotated the first de novo transcriptome of the head tissues of Salamandra salamandra from individuals, representative of distinct behavioral profiles. We obtained 1,153,432,918 reads, which were successfully assembled and annotated. The high-quality of the assembly was confirmed by three assembly validators. The alignment of contigs against the de novo transcriptome led to a mapping percentage higher than 94%. The homology annotation with DIAMOND led to 153,048 (blastx) and 95,942 (blastp) shared contigs, annotated on NR, Swiss-Prot and TrEMBL. The domain and site protein prediction led to 9850 GO-annotated contigs. This de novo transcriptome represents reliable reference for comparative gene expression studies between alternative behavioral types, for comparative gene expression studies within Salamandra, and for whole transcriptome and proteome studies in amphibians.

Unveiling a hotspot of genetic diversity in southern Italy for the endangered Hermann's tortoise Testudo hermanni.

BACKGROUND: Hotspots of intraspecific genetic diversity represent invaluable resources for species to cope with environmental changes, and their identification is increasingly recognized as a major goal of conservation ecology research. However, even for iconic and endangered species, conservation strategies are often planned without thorough information on the geographic patterns of genetic variation. Here, we investigated the spatial patterns of genetic variation of the endangered Hermann's tortoise Testudo hermanni in the Italian Peninsula by genotyping 174 individuals at 7 microsatellite loci, with the aim to contribute to planning effective conservation strategies. RESULTS: Ordination-based and Bayesian clustering analyses consistently identified three main genetic clusters, one spread in the central and northern part of the peninsula, and two restricted to southern Italy and Sicily, respectively. The highest levels of genetic diversity were found in populations of the southern cluster and, in particular, at the northern edges of its distribution (He > 0.6, Ar > 2.8 ), that correspond to areas of putative secondary contact and admixture between distinct lineages. Our results clearly identify a hotspot of genetic diversity for the Hermann's tortoise in southern Italy. CONCLUSION: We inferred the evolutionary history and the spatial patterns of genetic variation of the Hermann's tortoise in the Italian Peninsula. We identified three main genetic clusters along the peninsula and a hotspot of intraspecific diversity in southern Italy. Our results underline the urgent need for conservation actions to warrant the long-term persistence of viable tortoise populations in this area. Furthrmore, these data add further evidence to the role of southern Italy as a biodiversity hotspot for temperate fauna, claiming for higher consideration of this area in large scale conservation programs.

Brain de novo transcriptome assembly of a toad species showing polymorphic anti-predatory behavior.

Understanding the genomic underpinnings of antipredatory behaviors is a hot topic in eco-evolutionary research. Yellow-bellied toad of the genus Bombina are textbook examples of the deimatic display, a time-structured behavior aimed at startling predators. Here, we generated the first de novo brain transcriptome of the Apennine yellow-bellied toad Bombina pachypus, a species showing inter-individual variation in the deimatic display. Through Rna-Seq experiments on a set of individuals showing distinct behavioral phenotypes, we generated 316,329,573 reads, which were assembled and annotated. The high-quality assembly was confirmed by assembly validators and by aligning the contigs against the de novo transcriptome with a mapping percentage higher than 91.0%. The homology annotation with DIAMOND (blastx) led to 77,391 contigs annotated on Nr, Swiss Prot and TrEMBL, whereas the domain and site protein prediction made with InterProScan led to 4747 GO-annotated and 1025 KEGG-annotated contigs. The B. pachypus transcriptome described here will be a valuable resource for further studies on the genomic underpinnings of behavioral variation in amphibians.

From STRs to SNPs via ddRAD-seq: Geographic assignment of confiscated tortoises at reduced costs.

Assigning individuals to their source populations is crucial for conservation research, especially for endangered species threatened by illegal trade and translocations. Genetic assignment can be achieved with different types of molecular markers, but technical advantages and cost saving are recently promoting the shift from short tandem repeats (STRs) to single nucleotide polymorphisms (SNPs). Here, we designed, developed, and tested a small panel of SNPs for cost-effective geographic assignment of individuals with unknown origin of the endangered Mediterranean tortoise Testudo hermanni. We started by performing a ddRAD-seq experiment on 70 wild individuals of T. hermanni from 38 locations. Results obtained using 3182 SNPs are comparable to those previously obtained using STR markers in terms of genetic structure and power to identify the macro-area of origin. However, our SNPs revealed further insights into the substructure in Western populations, especially in Southern Italy. A small panel of highly informative SNPs was then selected and tested by genotyping 190 individuals using the KASP genotyping chemistry. All the samples from wild populations of known geographic origin were genetically re-assigned with high accuracy to the original population. This reduced SNPs panel represents an efficient molecular tool that enables individuals to be genotyped at low cost (less than €15 per sample) for geographical assignment and identification of hybrids. This information is crucial for the management in-situ of confiscated animals and their possible re-allocation in the wild. Our methodological pipeline can easily be extended to other species.

Mass of genes rather than master genes underlie the genomic architecture of amphibian speciation.

The genetic architecture of speciation, i.e., how intrinsic genomic incompatibilities promote reproductive isolation (RI) between diverging lineages, is one of the best-kept secrets of evolution. To directly assess whether incompatibilities arise in a limited set of large-effect speciation genes, or in a multitude of loci, we examined the geographic and genomic landscapes of introgression across the hybrid zones of 41 pairs of frog and toad lineages in the Western Palearctic region. As the divergence between lineages increases, phylogeographic transitions progressively become narrower, and larger parts of the genome resist introgression. This suggests that anuran speciation proceeds through a gradual accumulation of multiple barrier loci scattered across the genome, which ultimately deplete hybrid fitness by intrinsic postzygotic isolation, with behavioral isolation being achieved only at later stages. Moreover, these loci were disproportionately sex linked in one group (Hyla) but not in others (Rana and Bufotes), implying that large X-effects are not necessarily a rule of speciation with undifferentiated sex chromosomes. The highly polygenic nature of RI and the lack of hemizygous X/Z chromosomes could explain why the speciation clock ticks slower in amphibians compared to other vertebrates. The clock-like dynamics of speciation combined with the analytical focus on hybrid zones offer perspectives for more standardized practices of species delimitation.

Paso doble: A two-step Late Pleistocene range expansion in the Tyrrhenian tree frog Hyla sarda.

The Tyrrhenian tree frog, Hyla sarda, is an amphibian endemic to the Tyrrhenian islands (Western Mediterranean). Previous investigations of its Pleistocene evolutionary history suggested that it colonised the northern portion of its current range, through a spatial diffusion process from the Sardinia island, during the last glaciation. However, southern and northern portions of the species' range experienced markedly different climatic conditions during the Late Pleistocene, suggesting the possibility of an unusual two-step process of demographic expansion. Here, we use Bayesian phylogeographic approaches to locate the ancestral area in Sardinia and to characterise better the demographic component of this expansion event. These analyses located the ancestral area for H. sarda populations along the central-eastern coast of the Sardinia island, within an area previously shown to host suitable bioclimatic conditions for H. sarda populations throughout the Late Pleistocene. Historical demographic reconstructions clearly showed that a two-step process of demographic growth fits well the data, with northern populations expanding later than Sardinia populations. The harsher climatic conditions occurred in northern islands during the glacial epoch, as compared to Sardinia, likely delayed tree frog colonisation of northern territories, and the associated demographic growth.

Reconstructing hotspots of genetic diversity from glacial refugia and subsequent dispersal in Italian common toads (Bufo bufo).

Genetic diversity feeds the evolutionary process and allows populations to adapt to environmental changes. However, we still lack a thorough understanding of why hotspots of genetic diversity are so 'hot'. Here, we analysed the relative contribution of bioclimatic stability and genetic admixture between divergent lineages in shaping spatial patterns of genetic diversity in the common toad Bufo bufo along the Italian peninsula. We combined population genetic, phylogeographic and species distribution modelling (SDM) approaches to map ancestral areas, glacial refugia, and secondary contact zones. We consistently identified three phylogeographic lineages, distributed in northern, central and southern Italy. These lineages expanded from their ancestral areas and established secondary contact zones, before the last interglacial. SDM identified widespread glacial refugia in peninsular Italy, sometimes located under the present-day sea-level. Generalized linear models indicated genetic admixture as the only significant predictor of the levels of population genetic diversity. Our results show that glacial refugia contributed to preserving both levels and patterns of genetic diversity across glacial-interglacial cycles, but not to their formation, and highlight a general principle emerging in Mediterranean species: higher levels of genetic diversity mark populations with substantial contributions from multiple genetic lineages, irrespective of the location of glacial refugia.

Mapping the geographic origin of captive and confiscated Hermann's tortoises: A genetic toolkit for conservation and forensic analyses.

The illegal trade has been threatening tortoise populations worldwide for decades. Nowadays, however, DNA typing and forensic genetic approaches allow us to investigate the geographic origin of confiscated animals and to relocate them into the wild, providing that suitable molecular tools and reference data are available. Here we assess the suitability of a small panel of microsatellite markers to investigate patterns of illegal translocations and to assist forensic genetic applications in the endangered Mediterranean land tortoise Testudo hermanni hermanni. Specific allelic ladders were created for each locus and tested on several reference samples. We used the microsatellite panel to (i) increase our understanding of the population genetic structure in wild populations with new data from previously unsampled geographic areas (overall 461 wild individuals from 28 sampling sites); (ii) detect the presence of non-native individuals in wild populations; and (iii) identify the most likely geographic area of origin of 458 confiscated individuals hosted in Italian seizure and recovery centers. Our analysis initially identified six major genetic clusters corresponding to different geographic macro-areas along the Mediterranean range. Long-distance migrants among wild populations, due to translocations, were found and removed from the reference database. Assignment tests allowed us to allocate approximately 70 % of confiscated individuals of unknown origin to one of the six Mediterranean macro-areas. Most of the assigned tortoises belonged to the genetic cluster corresponding to the area where the respective captivity center was located. However, we also found evidence of long-distance origins of confiscated individuals, especially in centers along the Adriatic coast and facing the Balkan regions, a well-known source of illegally traded individuals. Our results clearly show that the microsatellite panel and the reference dataset can play a beneficial role in reintroduction and repatriation projects when confiscated individuals need to be re-assigned to their respective macro-area of origin before release, and can assist future forensic genetic applications in detecting the illegal trade and possession of Testudo hermanni individuals.

Detection and Phylogenetic Characterization of a Novel Herpesvirus in Sooty Terns Onychoprion fuscatus.

Since 2005, we have recorded annual episodes of alphaherpesvirus outbreaks in chicks of magnificent frigatebird Fregata magnificens on the Ile du Grand Connétable Nature Reserve in French Guiana. In 2009, we found sooty terns, Onychoprion fuscatus, that live sympatrically with frigatebirds, with visible clinical signs of a potential viral infection. To determine if the symptoms observed in sooty terns could be associated with an alphaherpesvirus previously identified in frigatebirds, we carried out molecular screening of samples collected from seven individuals. We identified and characterized a novel viral sequence from five birds. BLAST searches, pairwise nucleotide, and amino acid sequence comparisons, as well as phylogenetic analyses confirmed that the sequence belonged to the Herpesviridae family, of the Alphaherpesvirinae subfamily. We observed that it clustered with strains isolated from Podargidae (Caprimulgiformes), Columbiformes, and Falconiformes, but was distinct from the frigatebird herpesvirus. We have tentatively named it Onychoprion fuscatus alphaherpesvirus 1, (OfusAHV1). These two sequences, although found syntopic on the Ile du Grand Connétable, belong to two distinct alphaherpesvirus strains. Thus, the clinical symptoms showed by sooty terns do not likely result from a cross-species transmission event. Future work is needed to better characterize the virus and to investigate herpesvirus prevalence in healthy, free-ranging sooty terns, and to assess the impact of the virus on population viability.

Biogeographic history moulds population differentiation in ageing of oxidative status in an amphibian.

Regulation of oxidative status plays a substantial role in physiological ageing. However, we know little about age-related changes of oxidative status in wild animals, and even less about the role of population history in moulding ageing rates. We addressed these questions by means of a common garden experiment, using the Tyrrhenian tree frog Hyla sarda as the study species. This species underwent a range expansion from northern Sardinia (source) up to Corsica (newly founded) during the Late Pleistocene, and then the two populations became geographically isolated. We found that, at the beginning of the experiment, Sardinian and Corsican frogs had similar concentrations of all oxidative status markers analysed. One year later, Corsican frogs had higher oxidative stress and suffered higher mortality than Sardinian frogs. Our results suggest the intriguing scenario that population differentiation in rates of physiological ageing owing to oxidative stress might be an overlooked legacy of past biogeographic processes.

Longitudinal study of Amphibiocystidium sp. infection in a natural population of the Italian stream frog (Rana italica).

Mesomycetozoean-induced infections (order Dermocystida, genus Amphibiocystidium) in European and North American amphibians are causing alarm. To date, the pathogenicity of these parasites in field conditions has been poorly studied, and demographic consequences on amphibian populations have not been explored. In this study, an Amphibiocystidium sp. infection is reported in a natural population of the Italian stream frog (Rana italica) of Central Italy, over a 7-year period from 2008 to 2014. Light and electron microscope examinations, as well as partial 18S rDNA sequence analysis were used to characterize the parasite. Moreover, a capture-mark-recapture study was conducted to assess the frog demographics in response to infection. Negative effects of amphibiocystidiosis on individual survival and population fitness were absent throughout the sampling period, despite the high estimates of disease prevalence. This might have been due to resistance and/or tolerance strategies developed by the frogs in response to the persistence of Amphibiocystidium infection in this system. We hypothesized that in the examined R. italica population, amphibiocystidiosis is an ongoing endemic/epidemic infection. However, ecological and host-specific factors, interacting in a synergistic fashion, might be responsible for variations in the susceptibility to Amphibiocystidium infection of both conspecific populations and heterospecific individuals of R. italica.

A rapid rate of sex-chromosome turnover and non-random transitions in true frogs.

The canonical model of sex-chromosome evolution predicts that, as recombination is suppressed along sex chromosomes, gametologs will progressively differentiate, eventually becoming heteromorphic. However, there are numerous examples of homomorphic sex chromosomes across the tree of life. This homomorphy has been suggested to result from frequent sex-chromosome turnovers, yet we know little about which forces drive them. Here, we describe an extremely fast rate of turnover among 28 species of Ranidae. Transitions are not random, but converge on several chromosomes, potentially due to genes they harbour. Transitions also preserve the ancestral pattern of male heterogamety, in line with the 'hot-potato' model of sex-chromosome transitions, suggesting a key role for mutation-load accumulation in non-recombining genomic regions. The importance of mutation-load selection in frogs might result from the extreme heterochiasmy they exhibit, making frog sex chromosomes differentiate immediately from emergence and across their entire length.

Hybridization and extensive mitochondrial introgression among fire salamanders in peninsular Italy.

Discordance between mitochondrial and nuclear patterns of population genetic structure is providing key insights into the eco-evolutionary dynamics between and within species, and their assessment is highly relevant to biodiversity monitoring practices based on DNA barcoding approaches. Here, we investigate the population genetic structure of the fire salamander Salamandra salamandra in peninsular Italy. Both mitochondrial and nuclear markers clearly identified two main population groups. However, nuclear and mitochondrial zones of geographic transition between groups were located 600 km from one another. Recent population declines in central Italy partially erased the genetic imprints of past hybridization dynamics. However, the overall pattern of genetic variation, together with morphological and fossil data, suggest that a rampant mitochondrial introgression triggered the observed mitonuclear discordance, following a post-glacial secondary contact between lineages. Our results clearly show the major role played by reticulate evolution in shaping the structure of Salamandra salamandra populations and, together with similar findings in other regions of the species' range, contribute to identify the fire salamander as a particularly intriguing case to investigate the complexity of mechanisms triggering patterns of mitonuclear discordance in animals.